Exploring the Relationship Between Salmon Populations and Protected Areas

EDS 222 Final Project
MEDS
Author
Affiliation
Published

November 25, 2022

Background

Data

Data Wrangling

Geospatial

In order to calculate percent protected by watershed, I used spatial intersection to find the overlap between the protected areas in the CPAD database and each monitored watershed. I then divided the total area of overlap by the total watershed area and multiplied by 100 to find the % protected. I removed any protected areas established after 2009, as this was the first year of data I included in my analysis.

General

After finding percent protected for each watershed I added this back to the population count data. I selected only for adult population counts, as these were the most consistent metric of salmon population data collected. Some years adult population counts were estimated using different methods for the same population, I took the average when this occurred. Although the data set included data starting in 1981, data was collected inconsistently across populations, with more monitoring occurring as time increased. In order to prevent this from impacting my analysis, I chose to only analyze data within a set time frame which had consistent data collection and only include populations which had data that spanned this time frame. For steelhead the time frame I used was 2009-2018, and for coho I used 2010-2017. I also removed populations which never had any fish observed during the time period of my analysis. For ease of interpretation I set 2009 as year 0 for steelhead and 2010 as year 0 for coho.

Analysis

\(populationcount = B_0 + B_1year_t + + B_2percentprotected + B_3Year_t * percentProtected +E_i\)

Code
#stelhead regression
lm(populationcount~percent_protected + year + year:percent_protected, data = steelhead_all)
#coho regression
lm(populationcount~year + percent_protected + year:percent_protected, data = coho_final)

Results

Steelhead

term estimate std.error statistic p.value conf.low conf.high
(Intercept) 218.12 176.00 1.24 0.216 -128.06 564.31
percent_protected 4.23 3.35 1.26 0.207 -2.35 10.82
year 52.32 33.21 1.58 0.116 -13.01 117.65
percent_protected:year -0.99 0.62 -1.59 0.113 -2.21 0.24

Increase in percent protected has a positive effect on the number of steelhead. When year is 0 (2009) there will be on average 4 more steelhead for each one percentage increase in percent protected. Increase in percent protected decreases the slope of the relationship between year and population count. -0.99 is the difference in the effect of year on steelhead population count for every one increase in protected area. Time has a positive effect on the number of steelhead. When percent protected is zero on average there will be an increase of 52 fish per year.

Using a significance level of 0.05, \(B_1\) , \(B_2\), and \(B_3\) are not statistically significant (table 1) . My 95% confidence intervals include zero (table 1). I fail to reject the null hypothesis.

Coho

Increase in percent protected has a negative effect on the number of coho When year is 0 (2010) there will be on average 2 less coho for each one percentage increase in percent protected. Increase in percent protected decreases the slope of the relationship between year and population count. -0.46 is the expected difference in the effect of year on coho population count for every one increase in protected area. Time has a positive effect on the number of coho. When percent protected is zero on average there will be an increase of 41 fish per year.

Using a significance level of 0.05, \(B_1\) , \(B_2\), and \(B_3\) are not statistically significant (table 2) . My 95% confidence include 0. I fail to reject the null hypothesis.

term estimate std.error statistic p.value conf.low conf.high
(Intercept) 452.79 164.32 2.76 0.006 128.91 776.67
year 41.64 39.60 1.05 0.294 -36.41 119.69
percent_protected -1.72 3.08 -0.56 0.576 -7.80 4.35
year:percent_protected -0.46 0.74 -0.62 0.534 -1.92 1.00

Discussion and Limitations

Citations